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Abstract 

An efficient and powerful computational method is presented to synthesize a meta-material to specified electro- 
magnetic properties. Using the periodicity of meta-materials, the Finite Element Methodology (FEM) is developed to 
estimate the reflection and transmission through the meta-material structure for a normal plane wave incidence. For 
efficient computations of the reflection and transmission over a wide band frequency range through a meta-material 
a Finite Difference Time Domain (FDTD) approach is also developed. Using the Nicholson-Ross method and the 
Genetic Algorithms, a robust procedure to extract electromagnetic properties of meta-material from the knowledge 
of its reflection and transmission coefficients is described. Few numerical examples are also presented to validate the 
present approach. 


Introduction 

In the recent years there has been a great deal interest in studying the material (called meta-material [1-3]) hav- 
ing simultaneously negative permittivity and permeability. These meta-materials do not exist in nature and are artifi- 
cially formed by embedding properly designed metallic/non-metallic inclusions in binding mediums. In such a 
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material, the wave vector k and the Poynting vector P are anti-parallel causing reversal of some basic electromag- 
netic wave propagation characteristics such as Snell’s law and Doppler effects. The reversal of these basic propaga- 
tion properties has lead to some potential applications of these materials such as high resolution lenses, 
sub wavelength resonators and enhancement of antenna gains. One of the challenges faced by the researchers in this 
field is lack of a robust method to characterize these materials. The characterization of the meta-materials is carried 
out in two steps. First, the reflection and transmission properties of these materials are estimated /measured using a 
computational/experimental method. In the second step, a proper inversion process is used to extract effective elec- 
tromagnetic properties from the knowledge of reflection and transmission coefficients. Experimental procedures are 
too expensive and time consuming to explore all design space, and are is usually not preferred in early stages of 
design cycle. In the present paper, an efficient and powerful computational method is presented to characterize/esti- 
mate effective electromagnetic properties of meta-materials. In this method, first, the metallic inclusions of interest 
are modeled using a commercially available Computer Aided Design (CAD) package. Using the Finite Difference 
Time Domain method the reflection and transmission through these inclusions are then estimated. In the final step 
using an inversion procedure based on two approaches, Nicolson-Ross method [4] and Genetic Algorithm, effective 
properties of the selected inclusions are extracted. Few examples with known effective properties are presented to 
validate the present method. 


Analysis of Meta-Materials 

In this section, the development of computational methods based on the Finite Element Methodology and also 
Finite Difference Time Domain (FDTD) method are presented. Figure 1(a) shows a metallic inclusion in a split ring 
form and Figure 1(b) shows a dielectric medium embedded with a periodic arrangement of split ring shown in Figure 
1(a). It is assumed that the structure is infinite in the x-y plane and has finite width in the z-direction. It also assumed 
that the structure is illuminated by a plane wave with normal incidence as shown in Figure 1(b). 
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Figure 1(a) Split Ring 



Figure 1(b) Dielctric slab embedded with split rings 


Finite Element Method: To estimate the reflection and transmission coefficients of such a structure, a 
unit periodic cell as shown in Figure 1(b) is considered. Using the electric and magnetic symmetries, the 
unit cell can be terminated by electric and magnetic conducting boundaries at y = +a/2 and 
x = +b/ 2 , respectively. The electromagnetic field incident at the input port of the unit cell can be writ- 
ten as 

„ -ik n z * - ik n z 

Ein = yE 0 e , H\ n = -xE 0 /r\ 0 e (1) 


The electromagnetic field scattered in the regions Z < 0 and z > d can be expressed in terms of the Floquet modes. 
Using the FEM formulation, the weak form of wave equation to be solved for the electric field in the unit cell takes 
following form [5] 
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where E is the electric field in the unit cell volume V bounded by z — 0 and z = d planes, e is the vector 
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modal function for the Floquet modes, Y is the modal admittance for the p mode, T is the testing function, and 

S i , S', are the surface areas at the Z = 0 and z — d planes, respectively. From the solution of equation (2) the 

reflection and transmission coefficients for the structure shown in Figure 1(b) can be determined. 

FDTD Method: It is often required to estimate the reflection and transmission of EM waves through meta-material 
over multiple frequencies. The FEM approach presented above becomes computationally intensive when required to 
be repeated for each frequency of interest. The FDTD method provides an alternative and efficient way of estimating 
reflection and transmission through meta-materials over multiple frequencies. To perform these computations using 
the FDTD approach we used the commercially available CST Microwave Studio [6] computational tool. To validate 
the FEM code and the numerical results obtained using the CST software, a dielectric slab embedded with four split 
rings and capacitively loaded strip dipole as shown in Figure 2(a) is considered. The reflection and transmission coef- 
ficients obtained using both approaches are shown in Figure 2(b). The close agreement between the two numerical 



results confirms validity of the FEM formulation and the also the modeling procedure followed while using the CST 
Microwave Studio. 

Using these techniques the reflection and transmission characteristics of metallic inclusions of various shapes 
and sizes are computed and results will be shared with the audience at the time of presentation. 



Figure 2(a) geometry of split ring and 
capacitive loaded strips 
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Figure2(b) Reflection, transmission characteristics of 
split ring/strip geometry 


Extraction of Electrical Properties: 

Following the Nicholson-Ross [4] method to estimate electric and magnetic properties of meta-materials, it is 
assumed that the dielectric slab with the metallic inclusions is replaced by a homogeneous material slab having effec- 


tive permittivity ^ 8 
can be written as 
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and permeability ( pO J . The reflection and transmission coefficients of this effective slab 
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properties of an equivalent slab can now be defined as determination of 8 r , |T ; . using (3) from the knowledge of 
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The steps involved in extraction process are: 
(a) determine (r, t) using (4), 


(b) express 
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(c) with proper choice of n, parameters are given by = (( 1 + r)/( 1 - r))k , Z r = (( 1 - r)/( 1 + r))k r 

Due to ambiguity in selecting proper value of n, earlier researchers [7] assuming d « A. had simplified the calcula- 
tions in the step (b). However, the assumption d « A, limits application to very thin slab. In this work we have devel- 




oped a procedure to track correct value of n without making the assumption. The proper values of effective electric 
properties are also estimated using the Genetic Algorithm. 

Numerical Results 

In this section numerical data for the effective parameters of a slab embedded with split rings arranged in a peri- 
odic fashion, shown in Figure 3(a), are presented. First, the unit cell including a single layer of split ring is modeled 
using a CAD package. The transmission and reflection coefficients as a function of frequency are then computed 
using the CST Microwave Studio software. Using the extraction procedure described above the effective parameters 
are estimated and are shown in Figure 3(b). Using the extracted parameters the reflection and transmission coeffi- 
cients are calculated and compared with the numerical data obtained using the CST Microwave Studio. 


Thin Metallic Split Rings 
Arranged In Periodic 
Fashion 


Composite Slab (Infinite in X-Y) 
Loaded with Metallic Split 
Rings 



Figure 3(a) Geometry of a slab loaded 
with split rings 



Figure 3(b) Effective parameters of slab loaded 
with split rings. 


Conclusions 

A computational electromagnetic method using the FEM and FDTD approaches have been presented to estimate 
the reflection and transmission through a meta-material slab. Nicholson-Ross technique[4] and also Genetic Algo- 
rithm procedure in conjunction with the knowledge of reflection and transmission coefficients have been used to 
extract effective permittivity and permeability of a meta-material slab. The reflection and transmission coefficient of 
actual meta-material slab have compared well with the numerical data obtained using the extracted effective permit- 
tivity and permeability. Due to page limitation the numerical results for a metallic inclusion in the shape of split ring 
only has been shown here. The metallic inclusions of other shapes will be discussed at the time of presentation. 
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